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I. INTRODUCTION 

What are typical features of a quantum system whose 
classical limit is chaotic? This is a central question of 
so-called "quantum chaos" [2 [21 El SI |5], a branch of 
quantum physics that has been attracting a considerable 
interest since 1970's. Apparently, quantum systems show 
no trajectories, hence no Lyapunov exponents, Poincare 
sections or other signatures constitutional for the dis- 
tinction of chaos on the classical level. Instead, some 
genuinely quantum attributes of the system seem to ab- 
sorb the information on the regular or chaotic charac- 
ter of the classical dynamics. The best known examples 
are correlation properties of the spectra of energy levels. 
As a surprise, quantum systems with a chaotic classical 
counterpart show highly correlated quantum spectra, de- 
scribed within the theory of Gaussian matrix ensembles 
[6], while the spectra of systems that are classically reg- 
ular look more or less random. 

In recent years, alternative signatures of quantum 
chaos have been proposed, like the morphology of wave 
functions [2 121 E] , fluctuations of the scattering matrix 
[5], sensitivity to perturbations [7] etc. The research of 
these issues is by far not completed. Note that the ab- 
sence of an exact deflnition of chaos on the quantum level 
led to a proposal to use the term "quantum chaology" in- 
stead of quantum chaos [8j. 

The relation of spectral properties of chaotic quantum 
systems to those of Gaussian matrices was proposed by 
Bohigas, Giannoni, and Schmit in 1984 [9 . Since then, 
the conjecture has been tested in numerous concrete sys- 
tems and supported by several involved theoretical anal- 
yses. Recently, correlation properties of quant al spectra 
were rephrased into the language of stochastic time se- 
ries with ~ 1//" type of noise, the chaotic case being 
identifled with a = l [10 . 

In spite of this progress, some problems concerning the 
relation of the level statistics to classical chaos remain 
open. The following two questions, in particular, helped 
to guide the work presented in this article: First, if the 
classical dynamics exhibits abrupt transitions between 
dominantly regular and dominantly chaotic types of mo- 
tions with varying energy, to what extent does the level 
statistics within a single spectrum follow these changes? 



Second, since the quantization is not a unique procedure, 
does Bohigas' conjecture hold in all quantum realizations 
of the given classical system? 

The model we use to probe the above questions is 
the geometric model of nuclear collective motions [11]. 
Classical dynamics generated by this model was recently 
shown to exhibit an immense variability of the dynami- 
cal modes [T21 US] . The rise of ordered modes from the 
chaotic ones and their breakdown are phenomena ob- 
served at numerous places in the plane of energy ver- 
sus control parameter — and in a majority of cases they 
are not properly understood yet. Very similar properties 
have been assigned also to the interacting boson model 
[m [151 im 5 which is closely related to the geometric col- 
lective model. A detailed analysis of dynamical features 
of these models is an interesting and important subject 
in the context of nuclear physics. However, we consider 
these models to be very well suited also for more gen- 
erally oriented studies, such as those seeking answers to 
the above questions. 

The classical analysis of chaos in the geometric collec- 
tive model has been presented in Refs. [T21 [13]. In the 
present paper, we focus on the analysis of quantum prop- 
erties. We restrict ourselves to a subset of quantum levels 
with zero angular momentum, which makes the conflg- 
uration space effectively two-dimensional, in contrast to 
flve-dimensional space corresponding to general motions. 
It is then shown that the classical version of the model for 
zero rotations can be quantized in two physically mean- 
ingful ways. We solve the eigenvalue problem in both 
cases and compare the level statistics obtained, looking 
particularly into the regions where transitions between 
regular and chaotic dynamics take place. The possibility 
of changing the value of a classicality (Planck) constant 
enables us to populate the spectrum with variable density 
of quantum states, which is used for a global inspection of 
large energy domains and zooming in some finer details. 

The previous paragraph outlined the content of the 
present part of this paper, further referred to as Part 1. 
In the forthcoming part. Part II [17], the method in- 
vented by Peres [TH^ will be applied to the geometric 
model. The method, which exploits specific information 
on the structure of individual eigenstates, makes it possi- 
ble to draw the spectrum of a quantum system in a way 
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that allows one to visually allocate regular and chaotic 
domains. A close relation to classical dynamics, in par- 
ticular an analogy with the graphical method based on 
Poincare sections, will be demonstrated. 

The plan of the present part of the paper is as follows: 
In Sec, [m (and in Appendix [A]) , we introduce the model 
with its alternative quantizations and discuss some tech- 
nical issues related to its numerical solution. The method 
used to evaluate the spectral statistics is described in 
Sec. |III| In Sec. [IV| the results of the statistical analy- 
sis are presented and compared with the corresponding 
classical measures of chaos. Conclusions are contained in 
Seem 



II. THE MODEL 
A. 5D and 2D Hamiltonians 

In this section we introduce the Hamiltonian of the 
Geometric Collective Model (GCM) for zero rotations, 
J = 0, and provide two different ways of its quantization, 
which are referred to as five-dimensional (5D) and two- 
dimensional (2D) cases. 

The kinetic and potential terms of the GCM Hamilto- 
nian H = T^V, 
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are built from generalized complex coordinates a = ^ 
(with ji = —2, . . . , +2) and the corresponding conjugate 

(2) 

momenta tt = tt/^ , which are both spherical tensors of 
rank 2. Parameter K in the kinetic part has the meaning 
of mass, while {A, B^C} determine the form of the poten- 
tial. Note that [a x b]^^"^ stands for a coupling of general 
tensors a and b to angular momentum A. The Hamilto- 
nian is rotationally invariant since it contains only scalar 
couplings of coordinates and momenta. 

The model is usually approached via an expansion 
of the nuclear radius into spherical harmonics, R = 

Ro{l-\- X^A/i (^^t^^^^i^^^)^ with only the A = 2 terms taken 
into account; hence the name geometric model. However, 
the coordinates can also have different interpretations 
keeping only their quadrupole tensor character. The geo- 
metric Hamiltonian in the above form was introduced by 
Aage Bohr in 1952 [19]. A way of systematic construction 
of higher-order terms in both potential and kinetic parts 
of the Hamiltonian was presented by Gneuss, Mosel, and 
Greiner [20]. Several other types of potential have been 
considered in connection with shape transitions in nuclei. 
An overview of these potentials with relevant references 
and the corresponding quantum solutions can be found 
in Ref. t21j. 



Coordinates a satisfy the constraint a^^''* = (— 
and therefore contain five independent real variables. 
Two of these variables capture the intrinsic shape of the 
nucleus (with a quadrupole deformation) and the remain- 
ing three variables describe the orientation of the nucleus 
in the laboratory frame (they can be associated with the 
Euler angles transforming the lab frame to the intrinsic 
one). In the intrinsic frame, only the shape variables 
are relevant. These are connected with the two indepen- 
dent scalar combinations of a's in Eq. Q and are usually 
parametrized as follows: 
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/3cos7 = x, a/2 Re a||^2 = 7 = ^ (3) 
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{(^±1 = 0). This yields the potential in the form: 

V = AP^^ BP^ cos 37 + Cp^ . (4) 

Standard quantization procedure with tt^ = —iJ^-^- 
leads to the kinetic term [19] 
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where T^ot stands for a nontrivial rotational part of the 
kinetic energy, containing the derivatives with respect to 
Euler angles and coupling all five dynamical variables. 
Since we restrict ourselves only to nonrotating regimes, 
J = 0, we set Trot = 0. This can be seen as a projec- 
tion of the full 5D coordinate system into an effectively 
two-dimensional space describing vibrational degrees of 
freedom P and 7. Note, however, that this projection 
differs from the 2D case by a modified definition of the 
scalar product, namely: 
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where ^f^(/3, 7), with z = 1, 2, stands for two wave func- 
tions and ^f{P^ 7) represents a measure.^ Appendix [A| 

gives some details, see in particular Eq. (A8). The nor- 
malization condition (^^^|^^^) = 1 leads to another 
class of solutions than would apply in case of a genuinely 
2D Hamiltonian. 

For J = 0, however, the truly 2D scheme represents 
an alternative way of quantization. It means that we put 
the system into the intrinsic frame first and only then 
carry out the quantization. The kinetic term obtained in 
this way reads as follows: 
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which is nothing but the standard 2D kinetic energy ex- 
pressed in polar coordinates (r, ^) = (/?, 7). The scalar 
product in this case is defined in the usual way, therefore 
the space of solutions coincides with L^(R^). Note that 
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FIG. 1: Energy levels of the GCM Hamiltonian for A = -1, B = 1.09, C = 1, and h^/K = 25 • 10~^ drawn inside the y = 
section of the potential well. Energy is given in relative units. Panels (a) and (b), respectively, show even and odd states in the 
2D quantization, panel (c) corresponds to the 5D quantization. Levels associated with the wave functions in Fig. |2] are marked 
with their respective ordinal numbers. In both 2D cases, levels of the harmonic- well approximation are drawn on the left. 



the Hamiltonian given by Eqs. ([7| and Q is a generaliza- 
tion of the widely studied Henon-Heiles model [22]. Let 
us stress that in the nuclear physics context only the 5D 
quantization is correct. 

Both forms ([5| and ^ have the same classical limit 
for J = 0. The corresponding Hamiltonians H^^ and 
H^^ , respectively, with the common potential Q, enable 
one to study the impact of the quantization method on 
spectral correlations. Although individual energy eigen- 
values obtained in both quantizations are different, we 
may assume — as implicit in Bohigas' conjecture — that 
the spectral statistics remains essentially the same (af- 
ter correctly separating levels with different conserved 
quantum numbers in both cases). The validity of this as- 
sumption will be discussed in Sec. |IV[ Note that despite 
the 2D and 5D quantum Hamiltonians carry a clear phys- 
ical meaning, they represent just two options out of an 
infinite number of quantization possibilities. 



B. Numerical solution 

Both versions of the J = GCM Hamiltonian are diag- 
onalized numerically, using the eigenbases of a 5D or 2D 
harmonic oscillator. The oscillator Hamiltonian reads as 
= T* + Ksc, where Ksc = ^osc/^^ (with Aosc being 
an arbitrary positive constant whose optimal choice will 
be discussed later) and T* stands for the 5D or 2D ki- 
netic operator ([5| or ([7|, respectively (we denote • =5D 



or 2D). A general GCM Hamiltonian i^* = T* + V is ex- 
pressed as H* = i^*g^ + V', where V = V^ — V^sc- It turns 
out that matrix elements of in both oscillator bases 
can be expressed analytically, which makes the process 
of numerical diagonalization very efficient. Details and 
explicit expressions can be found in Appendix [A| 

The original 5D solution of the GCM possesses several 
implicit symmetries, namely 



*^D(/3,7) = *'^(/?,-7), 
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with an arbitrary wave function in 5D. These rela- 
tions arise from the ambiguity of the system's orientation 
in the intrinsic frame. On the other hand, solutions of 
the 2D model do not a priori satisfy such symmetries. 
If the spectra associated with both quantizations are to 
be compared, conditions ([8| and (9]) need to be imposed 
externally also to the 2D case. This is done by select- 
ing a subset ^En^ (where E stands for even in variable 
7) of the 2D oscillator basis in which the diagonaliza- 
tion is carried out (see Appendix [A]) . In addition, if we 
relax condition ([9| and require only the equality of ab- 
solute values of the wave functions involved, we can take 
into account another independent class of 2D solutions, 
namely the wave functions odd in variable 7, i.e. satis- 
fying ^^^(/3,7) = -\E^^^(/3, -7). These are obtained by 
diagonalization in the subset ^q^^^ of the 2D oscillator 
states (Appendix [A|). 
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FIG. 2: Probability densities derived from the wave functions for selected levels from Fig. [T] Columns (a), (b), and (c) show 
the 2D even, 2D odd, and 5D cases, respectively. All distributions are constrained by the three-fold symmetry of the potential. 
The dashed lines demarcate the kinematically accessible area at given energy; cf. the dashed lines in Fig. ^ The states in two 
upper rows belong to the region where the quadratic- well approximation of the potential is valid. The third and fourth row, 
respectively, depict examples of regular and chaotic states with higher energies. 



The Hamiltonian matrix expressed in the truncated 
oscillator basis has a band form. The band width is ap- 
proximately equal to the maximal value of the principal 
quantum number in the selected subset of basis states. 
This makes the diagonalization feasible even at relatively 
high dimensions. The convergence of solutions is checked 
by a visual inspection of the distribution of eigenvector 
components in the oscillator basis (a bad convergence is 
signaled by missing tails of the computed distributions) 
and/or by trial calculations using variable size of the ba- 
sis (we test the stability of computed eigenvalues against 
an increase of the dimension). Our procedure guarantees 
that the precision 5E of individual eigenvalues satisfies 
the condition 5E <C A£^, where is an average spacing 
between levels in the selected part of the spectrum. This 
is needed for the determination of the nearest-neighbor 
spacing distribution. 

In order to make the net spectra of converged eigen- 
solutions as large as possible, we optimize the oscillator 
parameter Aqsc that determines a characteristic scale of 
the basis wave functions. The procedure is based on the 



determination of the Aqsc value for which the trace of the 
GCM Hamiltonian in the truncated basis is minimal (the 
optimal choice of Aqsc is however lower than this value, 
as empirically verified for the parameter ranges studied 
here). Taking all these issues into account, we have found 
that on a common personal computer one can employ up 
to 10^ basis states and obtain up to about 5 • 10^ well 
converging eigensolutions (exact numbers still depend on 
the choice of external parameters). 

An example showing all three classes of solutions — i.e., 
2D even, 2D odd, and 5D — is plotted in Fig.[l] Here, pa- 
rameters of Eq. Q were chosen such that the potential 
has a minimum at 7^ 0. Although the spectra in the 
three panels of Fig.[T]look different, very close similarities 
become apparent when comparing also the corresponding 
wave functions. This is done for selected levels in Fig. [2] 
At very low energies, the system is fully regular since the 
minimum of the potential Q can be approximated by 
a quadratic well. Wave functions belonging to this re- 
gion are seen in the first and second rows of Fig. (2] On 
the other hand, the region of mixed dynamics is exem- 
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plified by wave functions in the third and fourth rows, 
which correspond to regular and chaotic cases, respec- 
tively. Indeed whereas wave functions in the third row 
exhibit regular behavior (the wave function is localized 
within an area following some specific classical periodic 
orbits [5 ), wave functions in the fourth row show diverse 
structures and cover the whole accessible area. 

It should be noted that the 2D and 5D cases dif- 
fer in the differential element needed to calculate the 
probability distribution in the x x y plane. In the 
2D quantization we simply have \'^'^^ {x,y)\'^dx dy = 
|^2^(/3,7)P|J|c^/?d7, where J = f5 is the Jacobian of 
the transformation from {x^y) to 7). The first two 
columns of Fig. |2] show just the squared modulus of 
the respective wave functions. In the 5D case, how- 
ever, all matrix elements contain the measure Af from 
Eq. (|6|. In the rightmost column of Fig. |2] we show the 
squared wave function 7)^ multiplied by a factor 

J^/J = I sin 37 1 . The resulting 5D density therefore 
vanishes where sin 87 = 0, so it is visually similar to the 
2D odd case, in spite of the condition ([9]). This is also 
why the ground-state density in Fig. j^Jc) has two max- 
ima in each potential well, although the wave function 
itself has no node. 

Let us finally briefiy remark on the choice of param- 
eters in this work. It is closely connected with the 
scaling properties of the GCM Hamiltonian discussed in 
Refs.jlSl [23]. In the classical case [13], only one of the 
parameters {A^B^C} determines the scale-independent 
behavior of the system, while the others and K can be set 
to +1 or ±1 for A. In the quantum case [2T, the classical- 
ity parameter = /K (whose changes can be viewed 
either as changes of the Planck constant, or as changes 
of the mass) constitutes the second independent param- 
eter of the model which cannot, in general, be scaled to 
unity. This parameter determines the absolute density of 
states. 

In the following, we take B as the principal control 
parameter and choose n to locate a sufficient number of 
levels into the energy region of interest. The remaining 
parameters are fixed to (A, C) = (—1,+!), which in the 
nuclear context corresponds to nuclei with stable ground- 
state deformations. For 5 = 0, the system is completely 
integrable (since in this case the Hamiltonian does not de- 
pend on 7 and the xxy angular momentum is an obvious 
integral of motions). Therefore, the value of parameter B 
represents the strength of a nonintegrable perturbation. 

Let us note that all the above- introduced quantities 
and parameters are considered here dimensionless. The 
conversion to a concrete scale requires to choose the same 
unit for energy E and parameters A^ B, C, and k,. In 
nuclear context, the unit of {A, B^C} is set by the form 
of the potential (e.g. the depth of the minimum). The 
appropriate value of (connected with the effective mass 
parameter K) can then be determined by adjusting the 
number of states in a certain interval (e.g. below E = 0). 



III. SPECTRAL STATISTICS 

According to Bohigas' conjecture [9j, chaotic systems 
exhibit strong correlations between levels, which result in 
an effect known as "spectral rigidity". The short-range 
component of these correlations is most clearly mani- 
fested in the distribution of the nearest-neighbor spacings 
(NNS), i.e. gaps between neighboring levels in a trans- 
formed (so-called unfolded) spectrum. In fully chaotic 
systems, this distribution is amazingly well approximated 
by the Wigner distribution, while in systems with regular 
classical counterparts the NNS distribution tends to be 
Poissonian. 

A suitable quantity allowing one to interpolate between 
the two limiting cases is the Brody parameter co [24 . It 
is defined through the distribution 



P{s; cj) = (cc; + l)aujs'^ exp [-a^jS^^-^) , 



(10) 
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where s is the spacing between adjacent levels in the un- 
folded spectrum and a factor resulting from the re- 
quired conditions P{s] uo)ds = 1 (normalization) and 
sP{s]u)ds = 1 (unfolding). Eq. ( p!Q| ) interpolates 
between the Poisson {uj = 0) and Wigner(cc; = 1) distri- 
butions, hence a value uo G [0,1] obtained from a concrete 
spectrum tells us where between order and chaos the ac- 
tual system is. In spite of an artificial character of this 
interpolation, it has been argued that the Brody distri- 
bution is capable of fitting the data generated by realistic 
systems with mixed dynamics [25] , 

We use the following procedure: Eigenstates obtained 
by the Hamiltonian diagonalization in a sufficiently large 
basis are split into groups (bins), each of them counting 
1000 consecutive levels. The standard polynomial un- 
folding procedure [26 is applied in each group, removing 
locally a smooth part of the level density and normaliz- 
ing the average spacing to unity. To obtain the Brody 
parameter, we use the identity following from Eq. (10) 



lnln[l - /(s;w)]"^ = Ino;^ + (1 +a;)lns. 



(11) 



Here, I{s;uj) = P{s';uj)ds' can be estimated from a 
cumulative histogram of spacings in the unfolded spec- 
trum. A simple linear fit of the expression on the left- 
hand side of Eq. (11) in the logarithmic scale of variable 



s yields the Brody parameter uj for each bin of levels. 
An example is shown in Fig. |3] The bins subject to the 
above evaluations increase in energy and overlap with 
each other (the shift between successive bins was set to 
100 levels). The average energy of each bin is drawn on 
the abscissa in the resulting dependence of u on E (cf. 
Figs. l6}|5| below). 



Although the linear fit (11) can be easily implemented 



and demonstrates the validity of the Brody approxima- 
tion over a broad domain of spacings, it may increase 
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FIG. 4: Regular fraction /reg of the classical phase space as a 
function of the control parameter B and energy £' in a vicinity 
of E ^ 0. The degree of chaos is coded in shades of gray, with 
white (black) corresponding to complete order (chaos). One 
may notice an island with increased regularity near B ^ 0.6. 
The number of bins in the B direction is much larger than 
that in the E direction. 



FIG. 3: An example of the linear fit illustrating the deter- 
mination of the Brody parameter from an unfolded spectrum 
of levels within a single bin [it belongs to the dotted line in 
Fig. ^h) at E ^ 1.33]. The left-hand side of Eq. ([iT} is 



denoted as T. The NNS distribution and its fitted Brody ap- 
proximation (with uj = 0.62) are shown (by the full line) in 
the inset. Dashed and dotted lines correspond to Poisson and 
Wigner distributions, respectively. 



the relative weight of large values of s in the calcula- 
tion of uj [27|. We have performed a numerical simu- 
lation showing that the value of the Brody parameter 
determined in this way may be decreased by an absolute 
value Acjsyst ~ —0.08, while the statistical error result- 
ing from finiteness of the sample of levels is estimated 
as Acc^stat ^ ±0.07. These uncertainties should be taken 
into account when evaluating the dependence of co on 
energy, see Sec. |IV[ 

We tried to implement also the new method [10 based 
on the 1//" noise in spectral fluctuations, with a G [1, 2] 
corresponding to spectra in between fully chaotic (a = 1) 
and fully regular {a = 2) limits. The advantage of this 
method lies in its simple and elegant formulation (with 
no explicit reference to random matrix ensembles) and in 
the fact that it simultaneously captures both short- and 
long-range spectral correlations. 

In our case, however, the results were not quite satis- 
factory. The reason — an insufficient statistics — may be 
anticipated to be present in a majority of systems in 
which the competition between regular and chaotic mo- 
tions quickly varies with energy. Indeed, the work [10] 
demonstrates the power of the 1//" method on the Rob- 
nik billiard that (as all billiard or cavity systems) ex- 
hibits a constant, energy-independent ratio between reg- 
ular and chaotic phase-space volume. This allowed the 
authors to average over a huge number of successive sets 
of levels and to get very precise results. In contrast, prop- 
erties of individual GCM trajectories cannot be trivially 
scaled with energy. This feature, which in fact repre- 



sents an important motivation for the detailed analysis 
of the present system, results in a significant increase of 
the statistical error of the deduced (energy-dependent) 
exponents a, in some cases even exceeding 30%. 



IV. RESULTS 

In this section, we will compare the quantum measure 
of regularity (Brody parameter) obtained in the way de- 
scribed in Sec. |III| with the corresponding classical mea- 
sure. We will focus on the influence of the different quan- 
tization schemes and on the dependence of results on the 
classicality parameter hc. 

In our previous work [121 113 ^ classical measures of 
chaos in the geometric model were studied. We solved 
the classical equations of motions for a number A^tot 
of trajectories with fixed energy E and then, using the 
method based on so-called alignment indices [28], classi- 
fied each trajectory as either regular or chaotic [T3^. In 
order to quantify the overall degree of regularity at given 
energy, we calculated a regular fraction of the phase space 
/reg = ^reg/^tot7 whcrc A^reg represents the number of 
regular trajectories in the sample. The regular fraction 
/reg takes values from (fully chaotic dynamics) to 1 
(fully regular dynamics) and it can be compared with an 
adjunct (1 — a;) of the Brody parameter. 

In the classical case, the dependence of /reg on en- 
ergy and control parameter B exhibits very complex non- 
monotonous behavior, which is for energies around E = 
depicted in Fig. [4j The following features of this depen- 
dence are worth mentioning: First, the system is regular 
for small values of 5, since B = represents a fully inte- 
grable limit of the model (see the end of Sec.|ll|). Second, 
a well-pronounced island of regularity in a wide range of 
energies is observed at 5 ~ 0.6. As shown in Ref. [16], 
this region is connected with the so-called regular arc 
of the interacting boson model [HI [15]. Third, increased 
values of the regular fraction are observed for some values 
of B ai E ^ 0. The absolute energy E = corresponds 



7 



reg 







8=0.62 


\ 


(a) 








Classical regularity 




FIG. 5: A comparison between the classical regular fraction 
/reg (panel a) and an adjunct (1 —u) of the Brody parameter 
(panels b and c) for B — 0.62 (the main maximum of regular- 
ity at ^ = 0). Results for different values of k in the 2D even 
case are shown in panel (b), while results for different types 
of quantization for At = 25 • 10~^ are presented in panel (c). 
A reasonable agreement of classical and quantum measures 
is observed. The jittering of individual curves is caused by 
statistical errors. 



to a local maximum of the potential Q at /3 = 0. 

Figures |5}|7| depict the dependence of both measures 
/reg and (1 — cj) on energy. The classical measure is al- 
ways shown in the uppermost panel. The value of the 
Hamiltonian parameter B in the three figures was cho- 
sen diS B = 0.62, B = 1.09, and B = 0.24, respectively, 
with regard to the location of some extremes of regular- 
ity at £^ = [12 . The energy range shown in the figures 
represents the most interesting region, lying between the 
domains of full regularity at very low and very high ener- 
gies. For energies just above the global minimum of the 
potential the system is entirely regular due to the valid- 
ity of the harmonic-well approximation. With increasing 
energy, the regularity suddenly breaks down and con- 
tinues falling sharply until it nearly reaches zero. After 
this stage, it takes off again, somewhat surprisingly, and 



exhibits several well pronounced peaks of highly regular 
motions, which are separated by valleys of more chaotic 
dynamics. These structures can be seen in Figs. [5]-[7[ 
For sufficiently high energies, not shown in the present 
figures, the regularity starts growing steadily toward the 
fully regular limit, following roughly a logarithmic de- 
pendence. This is connected with the dominance of the 
(3^ term of the potential at high energies [121 112 • 

The lower panels of Figs. [5}|7|show the corresponding 
quantum measure, the adjunct of the Brody parameter. 
In accord with the discussion in Sec. |III[ we estimate the 
absolute errors of the (1 - co) curves as —0.07 and +0.11. 
The reason for different error sizes in up and down direc- 
tions is the above-discussed systematic error, which tends 
to underestimate the value of co. On the other hand, a 
numerical error of /reg was estimated as < 5% [13 . If the 
errors are taken into account, the fluctuations observed 
in the lower panels of all three figures are smoothened out 
and the overall correspondence between /reg and {1 — uo) 
becomes rather good. In particular, the observed max- 
ima and minima of both curves coincide. It turns out 
that the Brody parameter tends to slightly overestimate 
the regularity — this being so even if the above system- 
atic error is considered. Note that a similar behavior 
was observed for the exponent a in the 1//" noise anal- 
ysis [To]. Indeed, there is no reason to expect that /reg 
and (1 —u) behave in exactly the same way. We only ex- 
pect a qualitative agreement, and that is fully confirmed 
in the present calculation. 

A remark is needed concerning the minimum seen in 
Fig.jsjb) just below E = for lower values of Hi. It has no 
apparent counterpart in panel (a). Indeed, this minimum 
is only an artifact connected with the resonance of P and 
7 vibrational energies at B = 0.62 [16]. The appearance 
of nearly equidistant bunches of levels in the resonance 
region causes a serious deviation of the NNS distribution 
from the Brody form, which results in a nonrealistic value 
of the Brody parameter. The discrepancy is localized 
only in a relatively narrow interval and gets washed out as 
hz increases (the incriminated energy region is populated 
by a decreasing number of levels; a similar effect was 
discussed in Ref. [29]). 

The dependence of the quantum results on the classi- 
cality constant k, is shown in Figs. [5|b) and|6|b). We 
observe that the value of does not affect the energy 
dependence of the Brody parameter. Mutual deviations 
of the curves for various are bound inside the standard 
error interval. Instead, the curves for distinct tz differ 
in the width of the displayed energy range. It has the 
following reason: Since a decreasing value of the classi- 
cality parameter raises the density of the spectrum, the 
plots for smaller hz are more detailed but cannot reach 
higher energies because of computational limitations of 
the diagonalization procedure. In our case, 3-10^ re- 
liable energy levels were calculated for each value of hz 
(the dimension of the diagonalized matrix being about 
two times larger) and the curves are cut at th e ce ntroid 
energy of the uppermost bin of levels (see Sec. III). 
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6=1.09 




/r= 0.25-10" 



10 



20 



30 



FIG. 6: The same as in Fig.js] but for B — 1.09 (a minimum 
of regularity for E — 0). Only the 2D even quantization is 
analyzed for different values of k. 



B=0.24 



Classical regularity 




FIG. 7: Classical (a) and quantum (b) chaotic measures for 
B — 0.24 (a minimum of regularity for ^ = 0). Only re- 
sults with At = 25 • 10~^ are shown in panel (b) for various 
quantizations (fluctuations are caused by statistical errors). 



The dependence of the quantum measure of chaos on 
the method of quantization is shown in Figs, [sjc) and 
[7|^b). As was demonstrated in Fig. [l] the spectra ob- 
tained by the three types of quantizations (we consider 
the 2D quantization scheme further split into the odd and 
even cases, see Sec. II A) differ from each other distinctly. 



The question concerning the validity of Bohigas' conjec- 
ture in all quantizations has been raised above. The an- 
swer is evident from the present results: The Brody pa- 
rameter for all quantization schemes exhibits essentially 
the same dependence on energy, the mutual deviations 
being fully within the range of standard errors. There- 
fore, we can report that full agreement with Bohigas' con- 
jecture is confirmed in the present model independently 
of the chosen quantization scheme. 



V. CONCLUSIONS 

We have studied the correspondence between classical 
and quantum measures of chaos in the geometric collec- 
tive model adopted from nuclear physics. In spite of its 
conceptual simplicity, the model exhibits enormous com- 
plexity of solutions, with a very sensitive dependence of 
the observed behaviors on external parameters and en- 
ergy. These features qualify the model for being a subject 
of detailed analyses of the competition between regular 
and chaotic modes of motions on both classical and quan- 
tum levels. 

Although the dependence of chaotic measures on some 
external control parameters has been extensively studied 
in quantum billiards, see e.g. Refs. [30[ |31], the depen- 
dence on energy has been marginalized so far. This is 
partly because billiard systems do not permit this kind of 
analysis — their chaotic features are always uniform in en- 
ergy. Even the studies based on "soft" potential systems 
have been so far focused mostly on the cases with a rela- 
tively simple energy dependence of chaotic measures, see 
e.g. Refs. [22,^, 32 . The collective Hamiltonians used 
in Refs. [H [HI [El and also those studied in Refs. [T2| [T3] 
provide a very different perspective. In this sense, the 
present work can be considered as complementary to the 
studies based on two-dimensional billiard systems. We 
have demonstrated that a tight connection between clas- 
sical and quantum measures of chaos, embodied in the 
well-known Bohigas' conjecture [9 , remains valid even 
under the condition of a strong energy dependence. 

Another important conclusion of our work is the ob- 
served independence of quantum chaotic measures on the 
method used to quantize the system. Since the defini- 
tion of quantum chaos is based on the system's semi- 
classical limit, it would be very surprising to find the 
opposite — i.e. statistical properties of spectra depending 
on the quantization. However, the question deserves to 
be tested. The present work is probably only a first step 
in this direction. 

Finally, our results show that the Brody parameter, 
despite of being often deprecated, represents a reason- 
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ably sensitive measure of chaos in quantum system. We 
nevertheless do not stop at this stage. In the forthcoming 
part of this paper [17], features of the geometric model 
are analyzed with the aid of the method proposed by 
Peres [18 . This method enabled us to study the break- 
down and reoccurrence of ordered quantal spectra with 
running parameter 5 in a much more detailed way. 
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APPENDIX A: HAMILTONIAN MATRIX 
ELEMENTS 

Here we expose the 2D and 5D oscillator wave func- 
tions used for the diagonalization of both versions of the 
GCM Hamiltonian for J = and write down explicit 
expressions for relevant matrix elements. There exists a 
tight connection between the 2D and 5D bases, which 
shows up particularly in the form of matrix elements. 
Note that our derivation of matrix elements is based on 
Ref. [33 , but an alternative algebraic approach was re- 
cently presented also in Ref. [34 . 

We employ the 2D and 5D oscillator bases in radial 
coordinates P and 7. The oscillator Hamiltonian is given 
by ^osc = ^* + ^osc/5^ where • =2D or 5D and T* is 
the kinetic energy ([7|) or ([5|. The oscillator eigenstates 
in the 2D and 5D cases are denoted as \nm) and 
respectively, where n and v represent the corresponding 
radial quantum numbers. These states are also eigen- 
states of the angular momentum associated with the ro- 
tations varying angle 7, therefore they carry good quan- 
tum number m corresponding to the 0(2) invariant in the 
2D case, or fi corresponding to the 0(5) invariant with 
J = in the 5D case [35]. Matrix elements of a gen- 
eral GCM Hamiltonian in these bases can be calculated 
analytically. 

Starting with the 2D case, the basis wave functions 
read as 

(/3,7|nm) ^ = R'Zm'r^il) (Al) 

with the radial and angular parts given by 



^?L(7) 



^1^(7) 




2kn\ 



(n + 3m)! 
sin 3m7 



m 



0,1,.. 



cos 3m7 



(A2) 



if m = 
if m = 1, 2, . . . 

Here, LJ^ denotes the Laguerre polynomial and k = 
\/2A^^^ /h. The angular part is written for both odd 
(O) and even (E) cases. Note that we have applied the 



symmetry condition ([8|, selecting only the states with 
angular momentum quantum numbers equal to multiples 
of 3. The states are normalized as follows. 



X r ^Tm^{l)^TM dl = Sm'mSn'n , (A3) 

Jo 

with X standing for O or E. 

Following Ref. [33] (where however a different normal- 
ization of Laguerre polynomials is used) , we can calculate 
the relevant matrix elements. First, the oscillator Hamil- 
tonian itself trivially yields 



2D I 



hQ{2n + 3m + l)Sn'n^^^ 



n n^ra ra 7 



(A4) 



where = ^/2A^^^JK . To calculate matrix elements 
of a general Hamiltonian H^^ , we need to know matrix 
elements of the individual terms in V = {A — Aqsc)/^^ + 
BP^ cos 37 + C/3^. The radial parts of these elements can 
be read off from Tab. |l] and the angular contribution is 
given by 



{n\ m + 1| cos 37|nm) 



if m = 

if m = 1,2,... 



(A5) 



(for other combinations of basis states, the matrix ele- 
ments vanish). 

Let us turn now to the 5D case. Following the same 
procedure as above, we express the wave function 

(/3,7kM) ^ *Ep(/3,7) = (7) (A6) 

with radial and angular components 



2v\ 



r (1/ + 3/Li ■ 



-fct {k0 



$5D 



(7) = 



/2/Li + l 



(cos 37) 



(A7) 



satisfying the normalization 



f 

Jo 



27r 



*M?(7)*f (7) |sin 37I = S^,^S,,, . (A8) 



(A9) 



In analogy with ( |A4[ ) we have 

{iy',fi'\H',^,\ufi) = hn {2v + 3^+1) V.Vp 



The radial matrix elements can be simply obtained from 
Tab. [l] after substitution n ^ m ^ ji on the left-hand 
side, and n^i^, m^/i + l/2on the right-hand side. 
The angular part reads as 



{v ^ M + 1| cos 37|i^/i) 



v/(2/i + 1) (2/i + 3) 



(AlO) 
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TABLE I: Matrix elements of the radial potential terms in the 2D oscillator basis. Matrix elements for other combinations of 
the oscillator states are zero. 







= (2n + 3m + 1) 




nm) 


= -k~^ J(n + 1) (n + 3m + 1) 


{n,m-\- 


nm) 


= A;"^/ + 3m + 3) (n + 3m + 2) (n + 3m + 1) 


(n- l,m + l|/3^ 


nm) 


= -3/c"^/^ V^n (n + 3m + 2) (n + 3m + 1) 


(n- 2,m + l|/3^ 


nm) 


= 3A;"^/^ y^n {n - 1) (n + 3m + 1) 


(n-3,m + l|/3^ 


nm) 


= -fc-3/2^n(n-l)(n-2) 


{nm\/3^ 


nm) 


= /c"^ [n (n - 1) + (n + 3m + 1) (5n + 3m + 2)] 


{n-\-l,m\f3'^ 


nm) 


= -2/c"^ (2n + 3m + 2) ^/{n + 3m + 1) (n + 1) 


{n-\-2,m\f3^ 


nm) 


= k~'^y/{n + 3m + 2) (n + 3m + 1) (n + 2) (n + 1) 
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